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We propose a new mass matrix ansatz: At the grand unified (GU) scale, the standard model (SM) 
Yukawa coupling matrix elements are integer powers of the square root of the GU gauge coupling 
constant e = ^qgu, multiplied by order unity random complex numbers. It relates the hierarchy 
of the SM fermion masses and quark mixings to the gauge coupling constants, greatly reducing the 
SM parameters, and can give good fitting results of the SM fermion mass, quark mixing and CP 
violation parameters. This is a neat but very effective ansatz. 
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In the standard model (SM), the smallness of the quark 
and charged lepton masses and the quark mixing angles, 
and the hierarchy among them has been a problem for 
decades of years. Altogether with the CP violation phase, 
they constitute 13 of all the 19 free parameters in the SM. 
These 13 parameters may have a common origin at ex- 
tremely high energy scales and open us a narrow window 
to the underlying theory beyond the SM. The latest ex- 
perimental values of the SM fermion masses [l| can be 
translated to the running masses at Mz with the routine 
and method of Ref. [5] , wich is given in Table HI The 
three quark mixing angles and the CP violation phase 
can also be easily obtained from the latest experimental 
values of the four Wolfenstein parameters of the CKM 
matrix 
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In the SM, the neutrinos are definitely massless. 

Stitching proper mass matrices is an effective way 
to reduce the SM parameters, explore the fundamental 
physics beyond the SM and give meaningful predictions, 
and many different structures of mass matrices have been 



TABLE I: The SM fermion running masses at Mz- 
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proposed and researched intensively [B-Q- What's more, 
many relations among the 13 parameters are always in- 
corporated within the fermion mass matrices [8l-[l0| , such 
as the GU bottom-tauon mass unification to;, = to,- and 
the formula for the Cabibbo angle tan^ 9c ~ rnd/ms- 

If all the Yukawa matrix elements are of order unity at 
extreme high energy scales, e.g. at the GU scale Mq\j, 
it is impossible to generate the hierarchical SM fermion 
spectroscopy naturally, either with the SM or the su- 
persymmetric (SUSY) renormalization group equations 
(RGEs) % Thus the mass matrices must be constructed 
hierarchically at some high energy scale. 

There is an extremely natural small parameter at hand 
to construct the hierarchical mass matrices: acu — 
0.04 = 0.2^. It is more convenient to use e = ^/ogu — 0.2 
instead of acu, as 0.04 is too small. We can guess the 
smallness and hierarchy of the SM fermion masses and 
the three SM gauge coupling constants have common 
origin at Mq\j: They are related by acu = Coin- 
cidentally, e = ^/ogu — A, one parameter in Wolfenstein 
parameterization of the CKM matrix 0] . We yet stress 
that the mass matrices are constructed by e = ^/ogu 
but not by A, which, reflecting the hierarchy of the quark 
mixings, can be deduced from proper hierarchical quark 
matrices. 

For example, at Mquj 



(2) 



where nu,c,t, nuc,cu,ut,tu,ct,tc are all non- negative real 
numbers. The elements must be multiplied by phase 
factors to generate CP violation, it is a more natural 
choice to multiply order unity random complex numbers 
[ISinHIl- Specifically, we propose the arguments of 
the complex numbers be distributed evenly, but the log- 
arithms of the magnitudes of the complex numbers be 
of normal distribution, where the logarithmic operation 
guarantees that the magnitudes of the complex numbers 
are positive definite. If all the exponent numbers of the 
mass matrices are only non-negative real numbers, with 
no other constraints, they will be viewed as input param- 
eters. Then there will be more input parameters than the 
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TABLE II: Mgu and acu for different Msusy- 



MsuSY 


Mt 


1 TcV 


5 TeV 


Mgu (10'^ GeV) 


1.87 
24.71 


1.42 
25.87 


1.11 
26.89 



13 SM plienomenological parameters to fit, wfiicfi is op- 
posite to tlie initial desire to reduce the SM parameters. 
Alternatively, if they are all non-negative integers, then 
these integers can be viewed as factors from fundamen- 
tal theories, just as the factor 3 in the Georgi-Jarlskog 
mass relation 8] (if we consider the matrices are con- 
structed by powers of acu, then the exponent numbers 
are non- negative half- integers) . For example, they can be 
deduced from F-theory GUTs [ll|. It is rational to sup- 
pose further that Hue = Ucu, Uut = Uct = ntc, etc., 
leaving only 6 independent non-negative integers for each 
mass matrix. In this way, the fermion masses, quark mix- 
ings and CP violation, altogether 13 parameters, can be 
possibly deduced from proper matrix structures of the 
single parameter s = ^/ogv, and the GU gauge coupling 
constant acu naturally generates the three SM gauge 
coupling constants. The free parameters in the SM are 
greatly reduced. 

Choi [l^ has shown that GUTs with contorted mul- 
tiplets of fermions are also available, e.g. pairing of 
quarks and leptons from different generations such as 
{u,d) ® {Vfj_,p) and pairing of different left- and right- 
handed fermions such as {v^, e.)L® I^r- General and con- 
torted Yukawa matrices may lead to different values of 
observable quantities. For simplicity, we do not consider 
matrices with contorted structures. 

The goodness of fit (GoF) is defined to be [l^ 

^ \ Qi cxpt / 

where qi caic and qi expt are respectively quantities from 
theoretical calculations and experimental observations. 

Using the minimal supersymmetric standard model 
(MSSM) RGBs, the Yukawa matrices are evolved from 
Mgu to the SUSY scale Msusy- After the SUSY parti- 
cles decouple, the matrices are evolved further to lower 
energy scales with the SM RGEs. The one-loop SM and 
MSSM /3-functions needed can be found in Refs. pTtflst. 
At Mz, the SM fermion mass, quark mixing and CP vi- 
olation parameters are extracted to fit the experimental 
values. For simplicity, we ignore the possible effects on 
these SM parameters from the decoupling of the SUSY 
particles, the running of the tan /? — Vu/vd and the tran- 
sition between different renormalization schemes, and as- 
sume the SUSY particles all decouple at a common scale: 
Msusy 0, 111- We will use mainly Msusy = 1 TeV in 
the following calculation. However, Msusy has influence 
on the values of AfsusY and some extent, which 

is obvious in Table HIl 



To run the matrices below Msugy, we must first de- 
termine tan/3. If we suppose j/t ~ yf, ~ ~ 1 at Mquj 
tan /? must be large. To obtain tan /? conveniently, the 
three Yukawa matrices can be approximated 
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we will get at Mz 

ruu = 2.14 MeV, = 
rud = 3.21 MeV, rris = 
me = 0.951 MeV, = 



^0 o\ 

0, (4) 
Vo l) 



with order unity random complex numbers multiplying 
the (3, 3) elements. When the matrices are evolved from 
Mgu to Mz and the third generation masses are com- 
pared to the experimental values given in Table HI we ob- 
tain reliably tan/3 55, ignoring the effect from A/sugy- 
The value of tan /3 is consistent with the solution given 
in Ref. [13 • 

We find that when all the exponent integers of the 
three SM Yukawa coupling matrices take proper values 
at MgU) it is possible to obtain good fitting results of 
the nine SM fermion masses. For example, when 
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0.599 GeV, to* 
55.6 MeV, = 



-- 176.9 GeV, 
3.08 GeV, 



= 136 MeV, nir = 1.84 GeV(.6) 



Compared with the experimental values in Table HI using 
Eq. (I3|), we can obtain a small GoF = 0.68 for all the 
nine SM fermion masses. When the experimental errors 
are taken into consideration, the GoF is even smaller: 
= 0.53. 

Matrices with vanishing elements at A/gu are also pos- 
sible choice 043 ■ The exponent integers of the vanishing 
elements are infinity. For example, if we take 

/£» £A 

^/(AfGu) ^ £4 , (7) 



the up- type quark masses at Mz are 

m„ = 1.55 MeV, = 0.597 GeV, mt = 176.9 GeV 

with a very small GoF — 7.1 x 10"'^ for the three 
masses. Examples of down-type quark and charged lep- 
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ton matrices with vanishing elements are as follows: 

/ Q 



with mass eigenvalues at Mz 
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\0 £2 



rud = 2.21 MeV, m,. = 55.4 MeV, rub = 3.08 GeV 

8.8 X 10^3 for the three 



and a very small GoF 
down-type quark masses; 
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with mass eigenvalues at Mz 
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me = 0.337 MeV, 125 MeV, = 1.77 GeV 

and a small GoF = 0-18 for the three charged leptons 
masses. There are no constraints on the structures of the 
mass matrices; in principle we can construct any textures 
that lead to well fitted fermion masses. We can even 
construct matrices at Mqu , whose first and second family 
fermion masses are both completely from mixing with 
other generations. Here is an example: 



/o 



Its mass eigenvalues at Mz are 



£6 £5> 

£2 
£2 1 



(10) 



m„ = 1.38 MeV, = 0.821 GeV, mt = 176.3 GeV 

with a very small GoF — 0.073 for the three masses. 

It is worth noting that Msusy has considerable effect 
on the SM fermion spectroscopy. All the masses above 
are obtained with AfsuSY = 1 TeV. For the matrices in 
Eq. IS]), if we take Msusy = Mt, 

TUu = 2.62 MeV, = 0.673 GeV, mt = 173.8 GeV; 

(11) 

if Msusy = 5 TeV, 

rriu = 1.80 MeV, = 0.543 GeV, = 179.3 GeV. 

(12) 

Compared the up-type quark mass eigenvalues with dif- 
ferent Msusy, tbe effect on the first and second genera- 
tion masses is obvious. This is a general phenomenology 
for the quark and charged lepton masses, whatever struc- 
ture the fermion mass matrices are at Mqu- 

To obtain well fitted quark mixing and CP violation 
parameters, we must pare suitable up- and down- type 
quark matrices. Not all paring of good fitting up- and 
down-type quark matrices lead to well fitted quark mix- 
ings and CP violation, such as the first two matrices of 



Eq. ([5]) which would give much larger sin6'i3 = 0.0136 
than the experimental value sin6'i3 = 0.00359. When we 
choose 
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we get at Mz 

sin 6*12 = 0.227, sin^ss = 0.0189, sin6li3 = 
|sin(5| = 0.651 



(13) 



0.00402, 
(14) 



with a GoF = 0.76 for the four parameters. Quark 
matrices with vanishing elements are also available. For 

example. 
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lead to 



sin6ii2 = 0.219, sin 023 = 0.0 1 86, sin6li3 = 0.00257, 
|sin(5| = 0.649 

at Mz with a GoF x^ ~ 0.89 for the four parameters. All 
the mass eigenvalues of the quark matrices in Eq. (|13[) 
and Eq. ()15p are well consistent with experiments. 

Msusy also has possible effect on the quark mixing 
and CP violation parameters. For the quark matrices in 
Eq. (HSl), if Msusy = Mt, 



sin 6*12 = 0.231, sin 023 = 0.0211, sin6li3 = 0.00458, 
|sin(5| = 0.651 (16) 

with a GoF x^ = 0.65 for the four parameters; if 
Msusy = 5 TcV, 

sin6ii2 = 0.223, sin6l23 = 0.0 1 73, sin 6*13 = 0.00359, 
|sin(5| = 0.651 (17) 

with a GoF x^ — 0.90 for the four parameters. Compared 
the four parameters in Eq. ([16]), ([14]) and ([TT]) . the effect 
of A/susY on sin and sin 023 is small, the effect on 
sin 013 is comparatively large, but | smS\ is hardly subject 
to AfsusY- It is almost the same when the quark matrices 
are of other structures. The hierarchical structure of the 
CKM matrix, which is reflected by the powers of A in the 
Wolfenstein parameterization, can obtain naturally from 
proper quark matrices, as we stated above. 
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In conclusion, we have proposed a new ansatz for 
the SM fermion mass and mixing generation: At Mqu 
the Yukawa couphng matrices are constructed by non- 
negative integer powers of e = ^/ogv, i-e. by non- 
negative half- integer powers of acu: with matrix ele- 
ments multiplied by order unity random complex num- 
bers. By all the examples given above, we have demon- 



strated that this is a neat but extremely effective ansatz: 
It not only relates the hierarchy of SM fermion masses 
and quark mixings to the gauge coupling constants, re- 
ducing the number of the SM parameters to the utmost 
extent, but can give well fitted SM fermion mass, quark 
mixing and CP violation parameters. It could also apply 
to the neutrino mass and oscillation problem. 
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